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Abstract 

The occurrence of poles in the 2tt propagator in a hot pion medium is studied. 
The domain of non-trivial solutions to the corresponding bosonic gap equation is 
investigated as a function of the chemical potential and temperature of the gas. 
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When two heavy nuclei collide at very high energy, the matter that is created 
differs qualitatively from that traditionally studied in nuclear and elementary par- 
ticle physics. In the initial stages of the collision, copious production of gluons and 
quarks in a large volume leads to a rapid increase in the entropy, and the distinct 
possibility of a new phase of matter characterized by deconfined partonic degrees of 
freedom. One therefore hopes that relativistic heavy ion experiments at Brookhaven 
and CERN will provide some insight into the structure of the QCD vacuum, decon- 
finement, and chiral symmetry restoration. 

The hot transient system eventually evolves into a gas of hadrons at high en- 
ergy densities, whose interactions may well obscure many potential signals of the 
deconfined phase. To study this problem, extensive use has been made of sequential 
scattering models, such as hadronic cascades [p], @, 0- In principle, these trace the 
evolution of the system from hadronization to freeze-out, and can assign features 
in the data to the very early dense phase. For example, calculations using the 
Boltzmann equation with bose statistics to simulate heavy ion collisions at beam 
energies of 200 A ■ GeV |3], |J show that if the initial pion density is large enough, 
the low transverse momentum component of their spectra will be enhanced during 
the evolution of the pion gas ||. Implicit in these studies is a thermalization rate 
governed by free space tt-tt interactions. However, recent investigations || [7| of 
the effect of finite phase space occupation on intermediate states suggest that the 
t-matrix may well be suppressed within the environment of a relativistic heavy ion 
collision (RHIC) ||. This quantitatively affects the conclusions one may draw from 
such simulations, and points directly to the need to better understand hadronic 
properties and interactions in RHICs. 

Here we wish to reconsider it-it interactions in the presence of a dense and hot 
pion gas along the lines of our previous approach || . We shall address in some detail 
the question of pion pair formation in the medium. As we shall see, the in-medium 
2tt propagator exhibits a pole between a lower and an upper critical temperature, 
signaling a possible instability with respect to pion pair formation. 
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To simplify the solution of the 2tt propagator, we consider a separable model for 
the free space it-it interaction ||, in which the 2tt state couples to an intermediate 
a (£ = 0,1 = 0) or p (£=1,1=1) meson. The corresponding scattering matrix may 
be derived from the following effective Hamiltonian 

K = H — n (Nn — 2 N a ) 

= ^2(uj k - fi) b\b k + ( fl a - 2//) a\a a 

k a 

+ ~ £ {KK{°i + ° -a) {a\W\k x k 2 ) + h.cj , (1) 

kik2a 

where the b k (b\.) are the pion annihilation (creation) operators labeled by momen- 



tum k with energy uo k = yk 2 +m 2 . The a a (crl a ) are the corresponding operators for 
the intermediate meson labeled by a = (a, P) (—a = (a, —P)), where a = (£, I) are 



the spin-isospin quantum numbers. The state's energy is il a = \j P 2 + M 2 , where 
P = k\ + k 2 is the total momentum, and M a the mass of the exchanged meson. We 
chose a Yukawa coupling 

(a | W \k x k 2 ) = V2 {2n a 2uo kl 2uo k2 )-^ 2 V a (k) (2) 



V a (k) = 4tt g a uj k \j M a v a (k) (3) 

Va{k) = [i +w£ff ' (4) 

where k is the relative momentum of a pion in the cm. frame of the pion pair with 
total momentum P. The constants g a , M a and k a are chosen as in Ref. 0, and 
thus our Hamiltonian reproduces the it-it phase shifts in free space. 

The introduction of a chemical potential in the effective Hamiltonian may be 
justified as follows: There is considerable evidence that the pion distribution at 
freeze-out in RHICs is near local thermal equilibrium, but not in chemical equilib- 
rium. Since the it-it interaction is essentially elastic up to 1 GeV, in the later stages 
of the evolution the inelastic collision rate is insufficient to achieve full equilibrium. 
As the system continues to expand toward freeze-out, the dominant number con- 
serving processes still thermalize the system, but a non-zero chemical potential is 
built up &§. 



At a formal level in a field theoretical scheme it is not possible to define a total 
pion number operator from the canonical pion field. In the context of our semi- 
relativistic effective Hamiltonian the total pion number J2% ^ Afc * s ^ e isospin 
index) does not commute with H. It is, however, conserved on time scales larger than 
the life-time of the intermediate meson: the basic process is two pion annihilation 
to an intermediate a or p, which subsequently decays back into two pions. In this 
spirit, we introduce a chemical potential \x for the pion. The uncorrelated (causal) 
27r propagator is then given by (at P = 0) 

G ™ {E > k) = { e- 2(^-1*) + i V ~ E + 2(^-(.) + i V } coth ^ ' (5) 

Using the .fT-Hamiltonian Eq. (|I|), we may rederive the in-medium cr-meson prop- 
agator, which, after Fourier transformation, reads 

D-\E) = E 2 — (M a -2p) 2 

- 2 (M a -2fj) - J -jL-p \(0\W\k,-k)\ 2 Gj n (E,k). (6) 

The pole condition D~ l (E) = for the existence of bound 2n states (with P = 0) 
in the sigma channel (£ = 0, / = 0) follows from Eq. (|]) with (Q) and (0): 

* = 2 (M ff -2,) r**» % h) „.,,j { : k -»\,, co^-» 



gl v Jo E 2 - (M CT - 2p) 2 E 2 - A(tu k - fi) 2 2T 

= F»AE) ■ (7) 

In Fig. 1 we show the function F^^{E) for a fixed pion chemical potential of 
100 MeV, at five different temperatures (solid lines). These curves intersect with 
\jg a ~ 486 MeV (horizontal dashes) below energies of 2{m- K — //) (the bound state 
domain) if the temperature lies between a lower and an upper critical value of 
T « 80 MeV (£ polc = 2{m n - //)) and T c « 137 MeV {E pole = 0), respectively. 

Several observations can be made here: (i) One always obtains a pole if the 
temperature lies between T (/x) (pole at E — 2{m 1T — fi)) and T c (/i) (pole at E — 0). 
Thus, at fixed /x, no matter how weak the potential g a is (provided it is attractive 
in the neighborhood of the 2ix threshold), one always obtains a pole in a range of 
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sufficiently high temperatures (for fermions at a sufficiently low temperature). In 
practice, T and T c will exceed sensible values for pions as soon as /x drops below 
~ 100 MeV, since they are increasing functions of //; (ii) for fixed g a , the pole 
position shifts downward with increasing temperature (for fermions: pole position 
moves up with increasing temperature); and (iii) the pole position also shifts down 
with increasing g a (as for fermions). As a function of temperature, we therefore see 
a behavior for bosons opposite to that for fermions. 

The fact that increasing temperature reinforces the binding is somewhat coun- 
terintuitive, but it is an immediate consequence of the coth factor associated with 
bose statistics in Eq. (|?p (fermions: tanh), which enhances the medium density, 
effectively reinforcing the 2-body matrix element for increasing temperature. How- 
ever, the bound state energy drops to zero at a certain critical temperature T c {p). 
For T > T c , a real solution ceases to exist. As usual, such a zero eigenvalue signals a 
phase transition point. We shall see shortly that for T ~ T c the system has become 
so dense that pairs (corresponding to the pole) begin to obstruct each other and 
dissolve in the medium. 

In analogy with the Cooper pole in fermion scattering, a pole in the two boson 
propagator also implies a non-trivial solution to the corresponding "gap equation." 
It is straightforward to derive such an equation from the Hamiltonian Eq. (fl|) 

a B 2 2 [d?q v*(k) v a {q) A q E q 

Ak = 8 * g -/(2^ M a — Iji 2E q C ° th 2T • (8) 

where v a is given by Eq. (Q), and the boson quasiparticle energy is E\ = (u k — ji) 2 — 
Afc. The ansatz A& = Sv a (k) reduces Eq. (§) to a non-linear equation for the gap 
strength 5. 

In spite of the formal similarities of Eq. (||) with the corresponding fermionic gap 
equation, there are important differences: For bosons, the is subtracted in 
(fermions: added), and the temperature factor is a hyperbolic cotangent (fermions: 
tanh). We have seen previously that at sufficiently high temperature (> T (/i)) 
there exists a pole in the 27r propagator, i.e., finite T favors binding, contrary to the 
fermion case. As we shall see, a similar feature occurs for the gap. At fixed /i, above 

5 



a temperature Tq(/x) a finite gap appears in the spectrum, which disappears again 
beyond a critical temperature which will turn out to be T c (/i). At temperatures 
larger than T c , the medium density increases to the point where the pairs dissolve 
(i.e. 5 = 0). 

Let us now study the domain of existence of a real, finite value of A in the fi-T 
plane. Rewriting the gap equation flD as 

4 = 2 r dkk > m . 1 ^Egpgg 

9l h M CT -2/i ^(uj k -n) 2 -5 2 vl{k) 2T 

= GWr(<5) , (9) 

it is clear that for fixed \x and T, G^t{S) increases monotonically from 5 = to 
5 m ax = m n — fi. Further, for fixed /i, G^t also increases with temperature. This 
behavior is shown in Fig. 2 for /x = 100 MeV and various values of T (solid lines). 
The pair potential 5 is given by the intersection of the solid lines with the horizontal 
dashed line. With decreasing temperature, it varies from zero at T = T C to <5 max 
in.. 1 1 at T = Tq. In Fig. 3 we show the gap strength 5 = A^/v a (k) as a function 
of m n — fi. For a given temperature, 5 vanishes at some upper critical chemical 
potential, and increases to 5 max at a lower critical chemical potential. 

The critical temperature T c for which the gap vanishes is determined by lineariz- 
ing Eq. (H), i.e., setting 5 = in Eq. (|9|). The resulting equation is identical to the 
pole equation (0) at E = 0, and thus the gap and the bound state both disappear 
at the same critical temperature T c . Eq. (|7|) is equivalent to what is known in the 
fermion the particle-particle RPA. Thus we see that in the boson case, too, 

one has the well known property that a phase transition to a pair condensate is sig- 
naled by the appearance of a zero eigenvalue in the RPA eigenfrequencies (Thouless 
criterion) flTT| . 

For fixed /i, as T decreases, 5 increases (see Fig. 2), until it reaches a maximum 
value of m n — fi at T = Tq. At this point Ek=o = 0, and the gap in the excitation 
spectrum has ceased to exist. Any further lowering of the temperature would yield 
a purely imaginary quasiparticle energy. In order to understand what happens, we 
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investigate the pion momentum distribution 

" t = <^ = exp( £ Jr)-i + ^ coth l (10) 

where v\ = [(u k — /j l )/(2E k ) — 1/2]. Since E k=0 — > for T — > Tq, we see that 
n k= Q diverges at T = Tq. The quasiparticles start to condense in the k = state in 
very much the same way as in the unpaired state when real particles condense as 
the chemical potential approaches m n for decreasing temperature. The interesting 
point in the paired case is that the condensation of quasiparticles happens at /i<m 7r , 
since it is counterbalanced by the finite value of 5. 

The region between the solid lines in Fig. 4 summarizes the domain of non- 
vanishing 5 in the \i-T plane. We see that the temperature Tq - where the gap 
closes up the excitation spectrum to zero at k = 0, and quasiparticles begin to 
condense as singles - is very close to the T (dot-dashed line) where the pole in the 
2tt propagator ceases to exist. The entire /i-T domain of non-trivial solutions to 
the gap equation is therefore intimately linked to the region where a pole exists in 
the 27T propagator. 

In summary, we have shown that finite temperature induces real poles in the 
2n propagator, even for situations where there is no 2tt bound state in free space. 
The situation is analogous to the Cooper pole of fermion systems, and we therefore 
studied the corresponding bosonic "gap" equation. This equation has non-trivial 
solutions in a certain domain of the /i-T plane. Such a region always exists, even 
in the limit of infinitesimally weak attraction. This is different from the T = case 



discussed by Saint James and Nozieres [12], where a nontrivial solution to the gap 



equation only exists when there is a two boson bound state in free space. Our study 
has to be considered preliminary. The final aim will be to obtain an equation of state 
for a hot pion gas within a Bruckner-Hartree-Fock-Bogoliubov approach. One of 
the major questions in this regard is how to avoid a collapse, since it-it interactions 
are primarily attractive with no need for a hard core repulsion down to very small 
distances. Also, the subtle question of single boson versus pair condensation must 
be addressed (see Ref. |12[] and references therein). In view of the strongly off-shell 



nature of the problem, the tt-tt interaction should be chosen to be more realistic 
than is the case here. Work on these topics is in progress. 

We are grateful to G. Bertsch, P. Danielewicz and M. Prakash for discussions. 
This work was supported in part by the U.S. Department of Energy under Grant 
No. DE-FG02-93ER40713. 
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Figure Captions 

Fig.l The function F^t of Eq. (§) versus energy, for a fixed chemical potential 
of \x = 100 MeV and several temperatures (solid lines). Intersections of the solid 
lines with the horizontal dashed line correspond to poles in the in-medium 

a channel 2ir propagator. 

Fig. 2 The function G^t of Eq. @ versus the gap strength 5, for /i = 100 MeV and 
several temperatures (solid lines). Intersections of the solid lines with the horizontal 
dashed line (l/g£) correspond to solutions of the gap equation (|]). At the lower 
critical temperature Tq the quasiparticle energy vanishes (5 = <5 max ), while the gap 
vanishes at the upper limit T c . 

Fig. 3 The gap strength versus m n — \i at several temperatures. 
Fig. 4 The area between the solid lines corresponds to the domain in the [i—T 
plane for which the gap equation (§) has non-trivial solutions. The 2ir propagator 
has a pole in the region between T (dot-dashed line, -E po i e = 2(m n — //)) and T c 

{E po l e = 0). 
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